Phase transitions for random states and a semicircle law for the partial transpose 
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For a system of A*" identical particles in a random pure state, there is a threshold ko = ko{N) ~ N/ 5 
such that two subsystems of k particles each typically share entanglement if fc > fco, and typically 
do not share entanglement if fc < ko. By "random" we mean here "uniformly distributed on the 
sphere of the corresponding Hilbert space." The analogous phase transition for the positive partial 
transpose (PPT) property can be described even more precisely. For example, for A'' qubits the two 
subsystems of size k are typically in a PPT state if A: < fci :— N/A~ 1/2 and typically in a non-PPT 
state if fc > fci . Since, for a given state of the entire system, the induced state of a subsystem is 
given by the partial trace, the above facts can be rephrased as properties of random induced states. 
An important step in the analysis depends on identifying the asymptotic spectral density of the 
partial transposes of such random induced states, a result which is interesting in its own right. 
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INTRODUCTION 

If all that we know about a quantum system is its di- 
mension n (the number of levels) and that it is well iso- 
lated from the environment, a reasonable model - or at 
least a reasonable first guess ~ for the state of the system 
is a unit vector selected at random from the sphere of an 
n-dimensional complex Hilbert space T-L. If the system 
interacts with some part of the environment, represented 
by an ancilla space T-La, the quantum formalism suggests 
as a model the so-called (random) induced state, obtained 
after partial tracing, over Hq, a random pure state on the 
space H'^'Ha. The same description applies if we are pri- 
marily interested in a subsystem of an isolated system, the 
setup that is addressed in the abstract. 

The above is just one example of how a random 
paradigm arises naturally in the quantum context. In 
the last few years probabilistic considerations have be- 
come a very fruitful approach in quantum information 
theory, the highlights being the fundamental paper [l| 
by Hayden, Leung and Winter and, more recently, Hast- 
ings's proof that suitably chosen random channels provide 
a counterexample to the additivity conjecture for classical 
capacity of quantum channels 

Although random states have been considered for many 
years, their properties (e.g., are they typically entangled?) 
remained elusive. In this note we describe a reasonably 
general way to handle such questions. Of course, the in- 
duced state p being random, we can not expect to be able 
to tell what p is. However, we may be able to infer some 
properties of p if they are generic (that is, occur with 
probability close to 1). As it turns out, for many natural 
properties a phenomenon of phase transition takes place 
(at least when dimH is sufficiently large): the generic be- 



havior of p "flips" to the opposite one when s := dun'Ha 
changes from being a little smaller than certain threshold 
dimension sq to being larger than sq. 

For simplicity, we will focus on the random induced 
states mentioned at the beginning of the Introduction. 
This leads (see 0, Q) to a natural family of probability 
measures on Viji), the set of states on H, where s, the 
dimension of the ancilla space, is a parameter. For speci- 
ficity, consider V. = ® and let us concentrate on 
two properties: entanglement and positive partial trans- 
pose (PPT). This choice is based, first, on the importance 
of these concepts and, second, on the differences in their 
respective mathematical features, which allow to present 
the diverse techniques needed to handle the problems. 

Concerning the importance aspect, we note that de- 
tecting and exploiting entanglement - originally discov- 
ered in the 1930's [H - is a central problem in quantum in- 
formation and quantum computation at least since Shor's 
work on integer factoring. Next, the positive partial 
transpose is the simplest test for entanglement (Peres- 
Horodecki PPT criterion, see 0,3) and is at the center 
of the important distillability conjecture , a positive an- 
swer to which would give a physical/operational meaning 
to the PPT property. On the other hand, from the com- 
putational complexity point of view, verifying the PPT 
property is easy (just check whether the partial transpose 
of the state p is positive semi-definite), while deciding 
whether p is entangled is a computationally intractable 
(NP-hard) problem 

In the special case when n := diuiH equals s = dim "Ho, 
partial tracing over the ancilla space Jia leads to the uni- 
form distribution on I?(C") (i.e., uniform with respect to 
the Lebesgue measure, or Hilbert-Schmidt volume, de- 
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noted by vol). More generally, when s ^ n, the eor- 
responding probability measure ^Xn,s has a density with 
respect to the Lebesgue measure on I?(C") which has a 
simple form Q 



dvol 
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where Zn^s is a normalization factor. Note that ([T]) de- 
fines the measure p,n,s (in particular) for every real s ^ n, 
while the partial trace construction makes sense only for 
integer values of s. If s < n, the measure fin,s is con- 
centrated on the boundary of I?(C"), but still can be 
described analytically. Another way to implement these 
measures is to start from the complex Wishart-Laguerre 
matrices W„^s x with s degrees of freedom) [11|, a 
classical ensemble in statistics and mathematical physics, 
and to normalize them to have trace 1. 

In spite of the explicitness of the formula ([ij, it is 
not easy to find - even approximately, and even for 
s ~ n — d"^ - the probability that a random induced state 
has PPT or is entangled. This is because these traits 
are not encoded in a simple way in the spectral proper- 
ties of p. It was shown in 1^ ~ via methods of high- 
dimensional probability ~ that the proportion of states 
(measured in the sense of Pn,m i-e., the Legesgue mea- 
sure) that are un-entangled, or separable, is extremely 
small in large dimensions. This was extended to the case 
when s = dim Ha is slightly larger than n = in 13|, 
while, on the other hand, it was proved in [l[ that ran- 
dom induced states on ® are typically separable 
when s is proportional to n? — d'^. The paper [l2| also 
established that un-entangled states are extremely rare 
even among PPT states (again, when s = n = d^). How- 
ever, even such simple question as ^'Does the proportion 
of the PPT states among all states go to as the dimen- 
sion increases? originally asked in [l^ . has not been 
rigorously addressed prior to the work that we describe 
in this note. The results we summarize go a long way 
in filling the gap s in understanding of the phenomena in 
question (see (l5l.[l6| for details and references). We show 
that the threshold between entanglement and separabil- 
ity occurs when s is roughly of order n^/"^ ~ d'^ , and that 
the threshold between NPT (i.e., non-PPT) and PPT is 
when s ^ 4n = Ad^ . 

The heuristics behind the consequences stated in the 
abstract is now as follows. If we have a system of N 
particles (with D levels each) which is in a random pure 
state, and two subsystems of k particles each, then the 
"joint state" of the subsystems is modeled by a random 
induced state on ® with d = and s = D^^^fc^ 
In particular, the relation k ~ N/5, or N ~ 5k, corre- 
sponds exactly to s = d^. A similar argument applies to 
the PPT property. 



Another consequence of the results is that, for a large 
range of parameters, when the ancilla dimension s is 
roughly between and d^, a generic random state is 
both PPT and entangled. Such states are bound entan- 
gled, or undistillable [9| and, in spite of being entangled, 
are useless for purposes such as teleportation or super- 
dense coding (cf. [l7|). However, since for small systems 
(2 ® 2 and 2 ® 3, see [ll, H^) PPT is equivalent to sep- 
arability, one is tempted to think that bound entangled 
states are an anomaly, and that the PPT property re- 
mains a good proxy for separability in higher dimensions. 
Our results imply that this heuristic becomes mislead- 
ing for large systems and that PPT and separability are 
quantitatively very different properties. 

While, as we postulated, random induced states form 
the most natural family of probability measures on 
I?(C"), our methods are fairly robust and allow handling 
of other random models. For example, another popular 
way to construct random states is to consider mixtures 
of pure states. Our analysis applies to this model as well: 
if Vn^s is the distribution of -j lV'i)('0i|i where [tpi) 

are independent uniform pure states, then all the results 
presented for the measures pn,s remain valid mutatis mu- 
tandis for Vr, a. 



THE RESULTS 

We recapitulate the setting: n = dim'H, i}} \s a. (ran- 
dom) unit vector uniformly distributed on the sphere of 
Ti^C, and p = tic is a random state on H whose 

distribution is denoted by Pn,s- Further, we assume that 
n = d^ > 1 and H = C^^C^; states on Ti will be consid- 
ered entangled, PPT etc., with respect to this particular 
splitting. For definiteness, the partial transpose F will be 
the transposition in the second factor, i.e., defined (by its 
action on product states) via (ti (g) T2)'" = ti (g) tJ. 

The first result describes the phase transition between 
generic entanglement and generic separability. 



Theorem 1 [16[ There exist effectively computable con- 
stants C, c > and a threshold function sq ~ so{d) satis- 
fying 

cd^ Said) s^Cd^log^d, (2) 

such that if p is a random state on ® distributed 
according to the measure Pd^^s cind if e > 0, then 

(i) for s ^ (1 — e)sQ{d) we have 

P{p is separable) ^ 2cxp{—c{e)d^) and 

(ii) for s ^ (1 + s)so{d) we have 

P{p is entangled) ^ 2 exp(— c(e)s), 
where c(e) > depends only on e. 
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Let us mention that our methods extend also to the 
multipartite setting and to "unbalanced" systems such 
as (8) C^, di ^ c?2 - see [3 for precise statements. 

The idea behind the proof of Theorem[TJ based on tools 
from high-dimensional convexity, is quite general and can 
be used to estimate thresholds for other properties of 
random induced states (beyond separability), provided 
the set of states with given property is a convex subset 
JC C 'D{T-L) and has some minimal invariance properties. 
However, in the special case of PPT we have a more pre- 
cise result. 

Theorem 2 [3l Let p be a random state on ® 



tially transposed random state on C^°Cg)C^°, when a — \ 
and a = 4 (sample size 1 in each case). 



Set si{d) = Ad^ and let 



distributed according to yLfp^ 
e > 0. Then 

(i) for s ^ (1 — e)si{d) we have 

P(/9 is PPT) ^ 2exp(-c(e)d2) and 
ill) for s ^ (1 + e)si{d) we have 

V{p is non-PPT) ^ 2 cxp(-c(e)s), 
where c{£) > depends only on e. 

As noted in [l5|, it is likely that the sharp estimate 
in (i) is of order exp (— c(£)(i"'); this conjecture leads to 
interesting large deviation problems for matrices . 

The proof of Theorem[2] (except for the exponential es- 
timates on the probabilities, which require a unified ap- 
proach common to both Theorems) depends on methods 
of random matrix theory and, specifically, on the follow- 
ing result that identifies asymptotic spectral density of the 
partial transpose of random induced states, and which is 
of independent interest. 

If A is a Hermitian matrix, we will denote by Ainax(j4) 
and Aniin(^) the largest and the smallest eigenvalues of 
A. If a e R and cr > 0, the semicircular distribution 
Msc(a,(T2) is the probability measure with support [a — 
2a, a + 2a] and density (2^^)-! ^4(t2 - {x - a)^. We 
then have 

Theorem 3 fisf Given a > 0, let pd be a random mixed 
state on C^f^C^ distributed according to Pd'^^[ad'^j ■ Then, 
as d tends to +oo, the eigenvalue distributions of p^ ap- 
proaches the deterministic measure /isc(i.i/a) t^f-c- fol- 
lowing sense: for any interval I C R, the proportion of 
eigenvalues of p\ inside the rescaled interval -^I con- 
verges (in probability) towards Psc(i,i/a){I)- 

Moreover, we also have convergence of the extreme 
eigenvalues XmaxicP Pd) and Xminid^p^) to respectively 
1 -I- 2l^foL and 1 — 2/y/a, the endpoints of the support 
o/Msc(i,i/a)- 

It is easy to numerically "check" the conclusion of The- 
orem [3] (this was first noticed in [l^l). For example. Fig- 
ure [1] shows sample distributions of eigenvalues of a par- 
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FIG. 1. Histogram showing distribution of the eigenvalues of 
, where p is a random state on C^" ® C^" chosen according 
to the distribution /i2500.2500 {a = I, top) or ^i2500,ioooo (q = 
4, bottom). In both cases the median eigenvalue is about 

2^=4x10-^ 



Because of the link between random induced states and 
the Wishart ensemble Wn,s, Theorem[3]holds also for that 
ensemble (real or complex, although it is the complex set- 
ting that is most relevant to the quantum theory); in that 
case the rescaling factor d^ is not needed. We emphasize 
that this is rather unexpected since the asymptotic spec- 
tral density of the Wishart ensemble itself is given by the 
Marchenko-Pastur distribution 21 1 . 



THE MATHEMATICS BEHIND THE RESULTS 

Although Theorems [T] and [H have similar statements, 
the tools used in their proofs are very different, which 
parallels the differences in the computational complexity 
of PPT vs. that of entanglement. However, combining 
all the tools is often necessary to obtain the full strength 
of the results. 

We first describe the proof of Theorem [l] which is of 
geometric nature and where the concept of mean width 
plays a central role. To present it, let us introduce basic 
concepts associated to a convex body K C R™ containing 
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the origin in the interior (see |22| for more baekground). 
The gauge of K is the function j| • j|x defined for x E R™ 
by 

IIxIIa' = inf{t ^ : x G tK}. 
The polar (or dual) body of X is defined as 

A'° = {?/eR'" : (x,2/)s^l VxeK}. 

If u is a vector from the unit sphere S'™^^, the 
support function of K in the direction u is hniu) := 
maxx£K{x,u) = \\u\\k°- Note that hxiu) + /lif (— u) is 
the distance between the two hyperplanes tangent to K 
and normal to u. The mean width of K is then defined as 



w{K] 



hjilu) da{u) = 



>da{i 



where da is the normalized spherical measure on S*™"^. 

In our setting, the relevant convex body is = 5°, 
where S is the set of mixed separable states on "H = 
(g) C^. The ambient space R™ is the space of self- 
adjoint trace 1 operators on H (hence m = d^ — 1), where 
the maximally mixed state plays the role of the origin. 
Since K° = {S°)° = S (the bipolar theorem) and since 
separability of p is equivalent to \\p\\s ^ 1, the crucial 
question is whether w{S°) is smaller or larger than 1. An 
analysis of this question leads to the following value of 
the threshold function appearing in Theorem [1] 

s„{d) = w{S°f. 

Assertions (i) and (ii) can then be derived from concen- 
tration of measure, and the heart of the proof is showing 
([2]), especially the upper bound. 

Determining the threshold so((i) requires finding the 
typical value of the gauge associated to S, computing 
which ~ as we mentioned ~ is a hard problem. We take 
an indirect route and find the order of magnitude of the 
threshold using the machinery of high-dimensional geom- 
etry, especially the so-called MM*-estimate. 

The MM*-estimate (see 22, 2^) is a general theorem 
which relates the mean width of a convex body and the 
mean width of its polar. While the abstract formula- 
tion may require an affine change of coordinates, in the 
present situation, because of the symmetries of S (invari- 
ance under local unitary conjugations), we can deduce 
via simple representation theory the inequalities 

1 ^ w{S)w{S°) s: C\ogd, 

where C > is a universal constant. Since w{S) can 
be estimated by standard techniques of high-dimensional 
probability [l2| , this allows to establish the order of mag- 
nitude of w{S°) (and hence of so('^)) up to polylog factors. 



As indicated earlier, the same scheme yields estimates 
for the thresholds corresponding to other properties in- 
cluding the PPT, but it does not allow to recover the 
precise order Ad^ appearing in Theorem [2l However, 
the latter result (except for quantitative estimates on the 
probabilities, which require further work, again based on 
the concentration of measure) follows readily from The- 
orem [3l which describes very precisely the spectrum of 
the partial transpose of a random induced state: the 
PPT condition is equivalent to Aniin(Pd) ^ 0, which is 
generic if 1 — 2/y/a > 0; similarly, XminiPd) < is generic 
if 1 — 2/y/a < - hence a = ^ = 4 is the critical 
value. In turn, to show Theorem [3] we use the moment 
method, a standard technique from random matrix the- 
ory. The idea is to identify the asymptotic spectral den- 
sity by computing its moments. This leads to problems in 
asymptotic combinatorics: the moments of semicircular 
distributions are given by the Catalan numbers, corre- 
sponding to the dominant combinatorial terms, while the 
statements about convergence of extreme eigenvalues are 
proved by refining the calculations and carefully estimat- 
ing contributions of lower order combinatorial terms. 



CONCLUSIONS 

We established that random induced states on "H = 
(E) exhibit a phase transition phenomenon with re- 
spect to the dimension s of the ancilla space. We exempli- 
fied the phenomenon on two properties: positive partial 
transpose, for which the threshold value of s is Ad^, and 
entanglement, for which the threshold is d^ (up to a poly- 
log factor). This allows to determine whether two sub- 
systems of an isolated system typically share (or typically 
do not share) entanglement when knowing only the sizes 
of those subsystems, and similarly for the PPT property. 
In fact, we provide a "black box" approach which applies 
to many natural properties of quantum states. Our re- 
sults motivate further study of the geometry of sets of 
quantum states, and that of large deviation behavior of 
some random matrix ensembles related to quantum in- 
formation theory. 

We expect the probabilistic methods to continue to 
play a major role in quantum theory. Indeed, the latter 
field usually involves high-dimensional objetcs; for exam- 
ple, the quantum analogue of a byte (a state on (C^)*^* 
- a qubyte, one may say) "lives" in a space of dimension 
2^^ — 1 = 65535. While this makes numerical schemes 
mostly impractical (the well-known curse of dimension- 
ality), randomness is boosted by the presence of many 
free parameters (one may call this phenomenon the bless- 
ing of dimensionality). The current level of understand- 
ing of these aspects of the theory is arguably comparable 
to that of combinatorics in the 1950's, when the power 
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of the probabilistic method [24| began to be appreciated 
and, subsequently, the study of random graphs became 
an intensive area of research. 
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